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Abstract 

We consider a braneworld model with bulk-brane energy exchange. This allows 
for crossing of the w = —1 phantom divide line without introducing phantom 
energy with quantum instabilities. We use the latest Snia data included in the 
Gold06 dataset to provide an estimate of the preferred parameter values of this 
braneworld model. We use three fitting approaches which provide best fit parameter 
values and hint towards a bulk energy component that behaves like relativistic 
matter which is propagating in the bulk and is moving at a speed v along the 
fifth dimension, while the bulk-brane energy exchange component corresponds to 
negative pressure and signifies energy fiowing from the bulk into the brane. We 
find that the best fit effective equation of state parameter Weff marginally crosses 
the phantom divide line w = —1. Thus, we have demonstrated both the ability of 
this class of braneworld models to provide crossing of the phantom divide and also 
that cosmological data hint towards natural values for the model parameters. 

1 Introduction 

Recent observations [H [2] of distant type la supernova (SnIa) provide strong indication 
of a rather unconventional cosmic expansion profile for our universe at late times. It is 
suggested that we are currently in an accelerated expansion phase [HIS], contrary to what 
we would expect in a matter dominated universe, where gravity tends to slow down the 
expansion. This result has had dramatic consequences on cosmological models, leading 
to the assumption of the existence of dark energy, an unknown component of the cosmic 
energy content which is responsible for the observed acceleration. Many models have been 
proposed to explain the nature of dark energy. These include new forms of matter and 
energy, like cosmological constant, scalar fields (quintessence), phantom fields, hessence, 
k-essence etc., as well as theories of modified gravity, extra dimensions and brane models 
[1 13 El Ej (see P [91 [lOl [m [12] and [I3] for a list of possible dark energy models). 
Apart from confirming the presence of dark energy, SnIa data also provide valuable 
information about its characteristics, namely the equation of state it obeys. Detailed 



analysis of various available datasets implies that a redshift dependent equation of state 
parameter w{z) crossing the phantom divide line w = — 1 is consistent with the data 
[13] and may be favored by some datasets. If such crossing is verified, it will render a 
number of cosmological models as non-viable candidates for dark energy, as it happens 
with pure scalar field models of phantom or quintessence fields. Obtaining models which 
can account for both the late times accelerated expansion, as well as an equation of state 
parameter with w = —1 crossing becomes thus a high priority. 

In [15] a model was presented based on brane cosmology, where dark energy is a 
manifestation of extra dimensions and the existence of bulk matter and energy exchange 
between the brane and the bulk. The model yields accelerated expansion at late times, as 
well as an effective equation of state parameter Weff, which crosses the Wejf = —1 line. 
Modelling the bulk content as a slowly moving perfect fluid, we showed that the crossing 
can occur without invoking matter that violates the weak energy condition, provided 
that the dark radiation field enters with a negative sign contribution. This requirement 
can be bypassed if a more general ansatz is adopted. A similar but qualitatively different 
from a physical viewpoint situation occurs in models which consider a dark matter - dark 
energy interaction (see [16] for some examples). In [T7| we departed from the bulk fluid 
interpretation and assumed that the bulk pressure and energy exchange contributions 
are the dominant ones, while the dark radiation field is negligible in the range < 2; < 1. 
We thus derived constraints on the allowed range of parameters for the bulk pressure 
and energy exchange terms by demanding the model to conform to the general expected 
temporal behavior [TTl [T8] . Our goal in this paper is to perform a more detailed analysis 
of the allowed range of parameters by fitting the model directly to the Gold supernova 
la dataset and checking whether the best fit values indicate accelerated expansion and 
w = —1 crossing, instead of requiring such a behavior a priori. In doing so, we will also 
keep the dark radiation term and see whether there are solutions which yield a positive 
sign for it. 

We should note here that the back-reaction of cosmological perturbations on large 
scales can mimic the effects of background accelerated expansion [19], thus altering the 
results of any analysis that focuses only on the background evolution. We have chosen 
not to include the effects of perturbations in our analysis for two reasons: 

• The set of equations of cosmological perturbations are known to be non-closed on 
the brane [20], i-e. information about the complete bulk evolution and dynamics is 
necessary in order to obtain a solution. Since the model we are considering doesn't 
originate from an exact bulk solution, but rather from a reasonable ansatz for the 
bulk matter content, a perturbative treatment cannot be particularly helpful. 

• It turns out that at late times, when the Hubble radius is much greater than the 
bulk radius of curvature, the results of perturbations appear to be equivalent to 
those obtained when treating a conventional 4D cosmology ^21j . 
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2 General Framework 



We will briefly review the model presented in llTj and we will recast its equations 
in a form which is more convenient for the SNIa fit procedure. The model consists of 
a four-dimensional brane imbedded inside a five- dimensional bulk space. Both the bulk 
and the brane are allowed to carry some energy content, while there may also be an 
energy exchange term between the two. Einstein's equations in this setup are 

GmN = RmN — -^GaInR = 4j^^3 '^MN ■ (1) 

where M is the 5D Planck mass. The energy- momentum tensor Tmn is 

Tmn = T^MN + ^A/iv - Gmn^ - 9i,v (y^{y) Kl^N , (2) 

The first and second term are the bulk and brane content respectively, while the third 
comes from the bulk cosmological constant and the fourth is the brane tension term. The 
metric ansatz we use to study the cosmology of a Friedmann-Robertson- Walker brane is 



Gmn = a\y,tHj , (3) 




where 7,^ is the metric for the maximally symmetric 3-space on the brane. For a bulk 
energy-momentum tensor T^^^ we shall take the general matrix 

-Pb P5 \ / PBn" -ri'Pr, 

T^'^^l = { PbS^, , Tif^ = PsaH, 1 . (4) 



-^P, Pb J V -n^P5 P^fe^ 

Note the presence of the off-diagonal component T5 = P5 signifying the flow of energy 
towards (or from) the brane. Note also the anisotropic choice Pb 7^ Pb, in general. 
Similarly, for a brane energy-momentum tensor Tj^\ we shall take the matrix 



T'K = ¥ I "0 pS; ) . tZ = ¥10 P-"-", I . (5) 
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Following [151 l22l 1231 121], we can use the junction conditions at the brane to derive the 
effective Friedmann equation on the brane that determines the cosmological evolution in 
four dimensions. Then, using an appropriate ansatz for the bulk pressure and the energy 
exchange term and assuming that the brane energy density is low, so we can neglect 
quadratic terms in p, we can exactly solve the system of equations to find the Hubble 
parameter. The ansatz we impose for the bulk energy-momentum tensor components 
which enter the cosmological equations on the brane is 

Pb = Da" , P, = F a^ . (6) 
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and it can be physically justified once we model the bulk content as a relativistic fluid, 
slowly moving along the fifth dimension. The case ^ = v corresponds to a slowly moving 
fluid with equation of state parameter = — 1 — /i/3. Note that more than one "fluid" 
component may be present. In such a case, Pb = P5 = ^ 

Then, it is possible that Di ^ Di, while Fi <^ Fi {i ^ l)and therefore, at late times 
Pb ~ -Di a^^ and P5 ~ F2 (^) a^^, accounting for different exponents. 

Substituting these expressions, we arrive at the effective Friedmann equation 

f k Srr ^ ,_, 

[a) + ^ = T^^^^^^' 

where Gn = 3'-// An = 3a /A'k{2AM^)^ is the AD Newton's constant and the effective 
energy density pejf stands for 

C C/27 26 ^ 2{3w - 1)F ^ 

P-ff - ^3{TT^ + - + + (/i + 4) [3{1 + w)+^if ■ ^ ' 

where 7 = a{2AM^)~'^ and we have defined 6 = D/2AM^. We have assumed the equation 
of state p = wp for matter on the brane. The constant C is just an integration constant 
for the energy density on the brane p(a). The integration constant C multiplies the 
dark radiation term and is related to global properties of the the bulk. The third and the 
fourth term in the expression for pe// come from the bulk pressure Pb and the brane-bulk 
energy exchange term P5 respectively. The values of the parameters C, F, u and /i, 
related to the dark radiation term and the bulk contributions are crucial in determining 
the profile of cosmic expansion. 

In [17] we performed an analysis for the allowed range of values for the latter four 
parameters, requiring a cosmological evolution that conforms to the general profile sug- 
gested by observational data, based on the assumption that the dark radiation term is 
negligible at late times. Here we shall perform a more general analysis, keeping also the 
dark radiation contribution and performing a fit to the Gold dataset, thus obtaining best 
fit values for the five parameters. In order to simplify the procedure, we will rewrite the 
effective Friedmann equation in terms of energy density parameters, assuming also a fiat 
space (fc = 0). It should be noted that actually there is aninteresting degeneracy between 
spatial curvature and time- varying dark energy as pointed out in Refs [251 ES] . Therefore 
a small spatial curvature term could wipe out the mild trend for w=-l crossing indicated 
by the GoldOG data in the context of a fiat Universe. However, we have chosen not to 
allow for a possibility of non-zero spatial curvature for two reasons. 

• There is significant theoretical motivation coming from inflation that k=0 to a 
very high accuracy. The possibility of non-zero spatial curvature would practically 
invalidate the flatness problem resolution provided by inflation. 

• Allowing for of non-zero spatial curvature would introduce one more parameter to 
be determined by the Snia data. However, the number of parameters involved in 
our brane world model is already large and this implies a large error region (see 
Fig. 1). The introduction of one more parameter would further extend this region 
making the results of our analysis very hard to interpret. 
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Doing so we obtain the following expression 

^ = 0-3(1+-) ^ ^DRa-^ + fisa" + fi5a^ (9) 

where we have defined the energy density parameters for brane matter, dark radiation 
and bulk components as 

"^<^-'-'"- 24M34> + 4) - 

and 

o .1 = SttG^ 2 (3^ - 1) F 

The brane energy density parameter corresponds to the observed matter density. 
It has a value of approximately Q^n — 0.27 [27j. Since the brane content is matter 
dominated, we will assume a value of w = in the above expressions. The four density 
parameters are in fact related by the fiatness requirement, which yields 

f^s = 1 — Qm — ^DR — ■ (13) 

Thus, we can eliminate the term in favor of ^2^. We see that it depends on both the 
coefficient D, which characterizes the bulk pressure, as well as the power u. In order to 
make the dependence on u explicit, we could rewrite it as = 

Using the above parametrization, we are now in position to fit the model to the su- 
pernova la data. There are four parameters to fit, Qdr, ^b, and /x. Having determined 
the best fit values for these parameter, we can infer the temporal behavior of the deceler- 
ation parameter q, as well as the profile of the effective equation of state parameter Wgff. 
We refer to it as an effective parameter, since it isn't directly related to some particular 
form of matter, but is representative of the accumulated effects of both bulk matter and 
energy exchange, together with dark radiation due to the presence of extra dimensions. 
The parameter Weff is according to the prescription of [28] 

"^^^^ " ^ 3 rflna ' ^^^^ 
where 6H^ = H^/Hq — Q^nO'^^ accounts for all terms in the Friedmann equation not 
related to the brane matter. Substituting the above parametrization, we obtain the for- 
mula 



from which we can determine the variation of w^jf with time. As we see, all time depen- 
dence comes for the bulk terms. The deceleration parameter q is given by 

1 a _ + 2fi^^a-^ - (z/ + 2) fi^a'' - (/i + 2) l^ga/^ 

In the following treatment, we will re-express q and We// in terms of the redshift z, 
instead of the scale factor a. The two are related by a = Oq being the current scale 
factor for our four- dimensional universe. 
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3 Gold dataset fit 



We are going to use the updated Gold dataset (Gold06) compiled by Riess et al. 
consisting of a total of 182 Snia, located at distances ranging within the interval 0.024 < 
z < 1.755. This dataset allows a fairly accurate reconstruction of the expansion history 
of our universe at late times. 

The Gold dataset is compiled from various sources analyzed in a consistent and robust 
manner with reduced calibration errors arising from systematics. It contains 119 points 
from previously published data plus 16 points with 0.46 < z < 1.39 discovered recently 
by the Hubble Space Telescope (HST). It also incorporates 47 points (0.25 < z < 1) 
from the first year release of the SNLS dataset [3] out of a total of 73 distant SnIa. 
Some supernovae were excluded [1] due to highly uncertain color measurements or high 
extinction Ay > 0.5, i.e. the dimming of the SN magnitudes resulting from dust. A 
redshift cut was also imposed at cz < 7000km/ s or z < 0.0233, in order to avoid 
the influence of a possible local "Hubble Bubble", i.e. a region with a higher Hubble 
parameter value due to the existence of a large local void [29]. Thus, a high- confidence 
subsample was defined. 

The above observations provide the apparent magnitude m{z) of the supernovae 
at peak brightness after implementing the correction for galactic extinction, the K- 
correction and the light curve width-luminosity correction. The resulting apparent mag- 
nitude m{z) is related to the luminosity distance d^lz) through 

mth{z) = M{M, Ho) + blogio{DL{z)) , (17) 

where in a flat cosmological model 

DAz) = {l + z) r dz' ^° 



H{z']ai,...,an) ' 

is the Hubble free luminosity distance {Hodi/c). The parameters ai, a„ are theoretical 
model parameters and M is the magnitude zero point offset, depending on the absolute 
magnitude M and on the present Hubble parameter Hq as 

M = M + 5%io(-r7^) + 25 = 
Mpc 

= M - 5/oc/io/i + 42.38 . (19) 

The parameter M is the absolute magnitude and is assumed to be constant after the 
above mentioned corrections have been implemented in m{z). 

The data points of the Gold06 dataset, after the corrections have been implemented, 
are given in terms of the distance modulus as 

fJ^obsizi) = niobsizi) - M . (20) 

The theoretical model parameters are determined by minimizing the quantity 

X (ai,...,an) = 2^ 2 ,9 , (21) 
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where ^ and cx^ ^ are the errors due to flux uncertainties and peculiar velocity dispersion 
respectively. These errors are assumed to be gaussian and uncorrelated. The theoretical 
distance modulus is defined as 



fith{zi) = mth{zi) - M = 5logio{DL{z)) + /iq , 



(22) 



where 



/io = 42.38 — 5logioh , 



(23) 



fJ'thizi) depends also on the parameters ai, used in the parametrization of H{z) in 
equation (fTSi) . 

To find the best fit parameters of the parametrization iQ we used the prior fiom = 
0.27 and followed three approaches which we label as minimization, marginalization- 
minimization and marginalization- averaging with results which are consistent among the 
three approaches. These approaches can be described as follows: 

• Minimization: We minimized eq. (1211) with respect to all four parameters {^Idr, 
VLbi v and /i). The steps we followed for the minimization of ( l2Ti) are described in 
detail in Ref. [30]. The advantage of this method is that it treats all parameters on 
an equal basis and gives direct information on all four of them. On the other hand 
it introduces significant degeneracy, thus leading to large error bars for the best fit 
parameter values. 

• Marginalization-Minimization: In this approach we marginalize over VLor & 
[—1.5,5] and VLb € [—10,2] and minimize with respect to v and /i. To choose the 
range in which we will marginalize the two parameters, since we have no prior 
information for them, we minimized for several pairs of reasonable values for v 
and /i and found a wide enough range that contains all the best fit VLdr and VLb- 
Also, we have checked that our results are not affected if we choose a larger range 
than the one we used to marginalize. The marginalized is defined as: 



This is minimized with respect to u and jj,. We thus obtain best fit values for the 
exponents v and /i which are directly relevant to the physical properties of the 
two components Pb and P5. This approach has less degeneracy at the minimum 
compared to the previous direct minimization but it does not provide an estimate 
for VLjjb and VLb which are treated as nuisance parameters. Such an estimate is 
needed in order to construct the best fit effective equation of state parameter 
Weff{z), defined by eq. f|T^ . In order to obtain this estimate we fix u and to the 
best fit values mentioned above and minimize with respect to Qdr and Qb- 

• Marginalization- Averaging: The best fit values of Qdr and Qb obtained after 
marginalization as described above were also verified by considering the average 
values of these parameters instead of their values that minimize x^- We thus define 
the average values of the two parameters as: 




(24) 




J e-xV2 cIQdr dClB 



•B 



(25) 
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where par = {VIdr, ^b}- 

The best fit parameter values obtained with each one of the above methods are 
shown in Table I. The first approach cannot be used to draw a safe conclusion regarding 
the model under consideration due to the large error-bars of the best-fit parameters. 
On the other hand, the marginalization seems to hint towards u ~ —3 and a negative 
Qb, which corresponds to a positive bulk pressure Pb (see eqs. and ffTOl) ). Also, a 
positive corresponds to P5 < which means that energy is flowing from the bulk into 
the brane. This is an intuitively reasonable result which may hint towards a relativistic 
matter component which is propagating in the bulk and is moving at a speed v along 
the fifth dimension [151 [13 • The dark radiation term has been shown to be related to 
the generalized comoving mass Ai of the bulk fluid and its best-fit parameter Qdr is 
positive [31] . This a desirable feature of this model, since a negative value would indicate 
a negative A4. It also ensures a well-behaved cosmological evolution at early times, 
although the validity of the effective Friedmann equation in this regime is questionable. 
Finally, as we can see in Table I, the two approaches, 'Marginalization-Minimization' 
and 'Marginalization-Averaging' are in good agreement. 

The corresponding forms of the equation of state parameter Weff{z) and the decelera- 
tion parameter q{z) are shown in Figs 1 and 2 for the first two approaches ('Minimization' 
and 'Marginalization-Minimization'). Despite the large error-bars, the best-fit forms of 
Weff{z) and q{z) are in good agreement between the two methods. The plots correspond- 
ing to the third approach ('Marginalization-Averaging') are practically identical with the 
Figs lb and 2b of the 'Marginalization-Minimization' method (see also Table I). 

A noteworthy feature of Figs la and 2a is the presence of "sweet-spots" at z ~ 0.2 
and z ~ 0.55. The presence of "sweet-spots" at different redshifts for different parame- 
terizations, usually polynomial, is something which has been studied previously (see [32] ) 
and is a consequence of the ansatz used. However, upon marginalization the error region 
of Fig. lb is significantly more smooth than the one of Fig. la. 



Table 1: The best fit parameter values obtained with each one of the three methods. 
Note that the quantities that were determined with the Averaging procedure of eq. (12^ 
were not assigned an error due to the nature of the method. 



Method 


v 






^B 


^^5 


Minimization 


0.75 ±5.09 


-3.55 ±5.77 


0.29 ±2.45 


3.39 ±0.76 


-0.27 ±2.57 


Marg.-Min. 


-3.0 ± 1.1 


-0.8 ±0.3 


0.49 ±0.25 


-1.00 ±0.49 


1.24 ± 1.23 


Marg.-Av. 


-3.0 ± 1.1 


-0.8 ±0.3 


0.52 


-1.05 


1.26 
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Figure 1: The effective equation of state Weff for the 'Minimization' method (Fig. la) 
and the 'Marginahzation-Minimization' approach (Fig. lb). The large error region of the 
in Fig. la is due to the significant degeneracy of the model while the decreased error 
region in Fig. lb is due to the fact that the values of u and fi were assumed fixed during 
the minimization with respect to ^Idr and Qb- 



4 Conclusion-Outlook 

We have used the latest Snia data included in the Gold06 dataset [2] to provide an 
estimate of the preferred parameter values of a braneworld model with bulk-brane energy 
exchange. The advantage of this model is that it has a clear physical motivation based on 
fundamental physics. Even though the large number of model parameters that require 
fitting introduces large error bars, the combination of different methods used for the 
fitting has provided interesting hints for the preferred parameter values. 

Due to the significant degeneracy of the model used, the first fitting approach pro- 
vided best fit parameter values with large error-bars (see Table I) and does not allow 
for safe conclusions to be drawn. On the other hand, the marginalization approach hints 
towards a bulk energy component that behaves like relativistic matter which is confined 
in the bulk and is moving at a speed v along the fifth dimension while the bulk-brane 
energy exchange component obtained corresponds to negative pressure, which means 
that energy is flowing from the bulk into the brane. Despite of these natural values of 
parameters, the best fit effective equation of state parameter Weff marginally crosses the 
phantom divide line w = —l.\t also reproduces accelerated expansion at late times, with 
a transition from acceleration to deceleration around z k, 0.5, consistent with previous 
analysis [3l [1] . This demonstrates both the ability of this class of braneworld models to 
provide crossing of the phantom divide, while the cosmological data hint towards natural 
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0. 0.25 0.5 0.75 1. 1.25 1.5 1.75 0.25 0.5 0.75 1. 1.25 1.5 1.75 



Figure 2: The deceleration parameter q for the 'Minimization' method (Fig. 2a) and the 
'Marginahzation-Minimization' approach (Fig. 2b). The large error region of the in Fig. 
2a is due to the significant degeneracy of the model while the decreased error region 
in Fig. 2b is due to the fact that the values of z/ and /i were assumed fixed during the 
minimization with respect to ^dr and VLb- 

values for the model parameters. 
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